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Abstract 



Consider the initial-boundary value problem for a strictly hyperbolic, genuinely nonlinear, 
Temple class system of conservation laws 

_ [ U{t, a) = Ua{t), 

Ut + f{u)^^0, u(0,x) = u{x), < ^ (1) 



{t,b)^Ut{t), 

on the domain Q, = {{t, x) £ M.^ : t > 0, a < x < b}. We study the mixed problem (1) from 
the point of view of control theory, taking the initial data u fixed, and regarding the boundary 
data Ua, Ub as control functions that vary in prescribed sets Ua, Ub, of L°° boundary controls. 
In particular, we consider the family of configurations 

A(T) = {u{T,-) ; u is a sol. to (1), Ua eUa, Ub eUb] 

that can be attained by the system at a given time T > 0, and we give a description of the 
attainable set A{T) in terms of suitable Oleinik-type conditions. We also establish closure 
and compactness of the set A{T) in the topology. 
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1 Introduction 



Consider the initial-boundary value problem for a strictly hyperbolic, genuinely nonlinear, 
system of conservation laws in one space dimension 

ut + f{u)^^0, (1.1) 

u{0,x) = u{x), (1.2) 

u{t,a) = Ua{t), (1.3) 

u{t,b) = Ub{t), (1.4) 

on the strip fl = {{t,x) G ; t > 0, x £ [a, b]}. Here, u — u{t,x) £ R" is the vector of 
the conserved quantities, Ua, Ub are measurable, bounded boundary data, and the flux function 
f : U t-^ R" is a smooth vector field defined on some open set U C M", that belongs to a class 
of fields introduced by Temple (26[ for which rarefaction and Hugoniot curves coincide. We 
recall that, for problems of this type, classical solutions may develop discontinuities in finite time, 
no matter of the regularity of the initial and boundary data. Hence, it is natural to consider 
weak solutions in the sense of distributions. Moreover, since, in general, the Dirichlct conditions 



(L3)-(L4) cannot be fulfilled pointwise a.e. (see ||6|, |l^), different weaker formulations of the 
boundary condition have been considered in the literature (see 0, U] and references therein) . 
Here, following F. Dubois, P.G. LeFloch |l^, we will adopt a formulation of ( |1.3| )-( |l^ ) based on 
the definition of a time-dependent set of admissible boundary data, that is related to the notion 
of Riemann problem. 

In the present paper, having in mind applications of Temple systems to problems of oil reservoir 
simulation, multicomponent chromatography, as well as in traffic flow models, we study the effect 



of the boundary conditions (L3)-(L4) on the solution of (l.l)-( |l.2| ) from the point of view of 
control theory. Namely, following the same approach adopted in [H, Q for scalar conservation 
laws, we fix an initial data u G L°°([a, b]) and we consider the family of configurations 

A{T; Ua, Ub) = {u{T, •) ; w is a sol. to (1.1) - (Q), Ua G Ua, Ub G Ub] (1.5) 

that can be attained at a given time T > by solutions to (pTl|)-( pr^ ), with boundary data 
Ua, Ub that vary in prescribed sets Ua, Ub C L°°(]R+) of admissible boundary controls. In the 
case of scalar, convex conservation laws, it was proved in Q, by using the theory of generalized 
characteristics |l^, that the profiles w{x) which can be attained at a fixed time T > are only 
those for which the map x i— > ^ (M^) jg j^qj^ increasing. Under the assumption that f'{u) > 
for all u, and for solutions of the mixed problem ( |l.l| )-( |T^ ) on the region fl, this condition is 
equivalent to the Oleinik-type inequalities 

^"-(^^ ^ O^^^It^ for a.e. XGK5], (1.6) 

[D^w denoting the upper Dini derivative of w). For general nx n systems, a complete character- 
ization of the attainable set does not seem possible, due to the complexity of repeated wave-front 
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interactions. However, in the particular case of Temple systems, wave interactions can only change 
the speed of wave-fronts, without modifying their amplitudes, due to the special geometric fea- 
tures of such systems. Therefore, the only restriction to boundary controllability is the decay 
due to genuine nonlinearity. We then consider here a convex, compact set T G U, and provide a 
description of the attainable set 

A{T)=A{T;14°°,U°°), = L°°([0,T],r), 

in terms of certain Olcinik-type conditions. We also establish the compactness of A{T) in the 
topology. Applications to calculus of variations and problems of optimization (where the 
cost functional depends on the profile of the solution at a fixed time T) motivate the study of 
topological properties of A{T). 

The paper is organized as follows. Section ^ contains the basic definitions and the statement 
of the main results. We also provide in this section a review of the existence and well-posedness 
theory for the mixed problem (|l.l|)-(pr^), and a description of a front tracking algorithm that 
will be used throughout the paper. In Section ^ we establish some preliminary estimates, and 
a regularity result concerning the global structure of solutions to the mixed problem (|L1|)-(|L^ 
generated by a front tracking algorithm. The proof of the main results is contained in Section |4[ 



2 Preliminaries and statement of the main results 



2.1 Formulation of the problem 



Let f : U ^ M" be the flux function of the strictly hyperbolic system (1.1) defined on a 
neighborhood of the origin U C M". Denote by Ai(u) < ••• < A„(m) the eigenvalues of the 
Jacobian matrix Df{u), and let {ri(u), . . . ,r„(u)} be a basis of right eigenvectors of Df{u). By 
possibly considering a sufficiently small restriction of the domain U, we may assume that the 
following uniform strict hyperbolicity condition holds. 

(SHI) For every u,v ^ U, the characteristic speeds at these points satisfy 

Xi{u)<Xj{v), yi<i<j<n. (2.1) 



We also assume that there is a fixed set of characteristic lines entering the interior of the strip 
[a, b] X K"*" at the boundaries x = a, x = b, i.e. that, for some index p G {1, . . . , n}, there holds 

Xp{u) < < Xp+i{u), \/ ueU, (2.2) 

and we let A""'^, A™'*'' denote the minimum and maximum characteristic speed so that there holds 

< A'"'" < |A,(u)| < A'""'^, yu€U. (2.3) 
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Moreover, we assume that each i-th characteristic field r.i is genuinely nonlinear in the sense 
of Lax 1^ , and that system (^]|) is of Temple class according with the following. 

Definition 2.1 A system of conservation laws is of Temple class if there exists a system of 
coordinates w = (wi,...,w„) consisting of Riemann invariants, and such that the level sets 
{u £ U; Wi{u) — constant^ are hyperplanes (see /p^J. 

By possibly performing a translation of coordinates, it is not restrictive to assume that the 
Riemann invariants are chosen so that di Xi{w) > 0, i = for all w = w{u), u £ U. 

Throughout the paper, we will often write Wi{t,x) = Wi(u{t,x)^ to denote the i-th Riemann 



coordinate of a solution u = u{t, x) to (1.1). We recall that, for a Temple class system, Hugoniot 
curves and rarefaction curves coincide . Moreover, as observed in , thanks to the existence 
of Riemann coordinates one can show that the assumption SHI implies the invertibility of the 
map /:f/^ /([/). 

We next introduce a definition of weak solution to ( |1.1| )-( [L^ ) which includes an entropy 
admissibility condition of Olcinik type on the decay of positive waves, so to achieve uniqueness. 



The boundary conditions (1.3)-(1.4) are formulated in terms of the weak trace of the flux /(w) at 
the the boundaries x = a, x = and are related to the notion of Riemann problem in the same 
spirit of [|l8j. To this purpose, letting u{t,x) — — x/t\ U]^,Ufj), u^, £ J7, denote the 
self-similar solution of the Riemann problem for (^]^) with initial data 

{ul if a; < 0, 
uh if a; > 0, 

for any given boundary state u £ U, we define the set of admissible states at the boundaries 

Va(u) = {W{0+; u,ur) ; u„, £ U}, 

(2.4) 

Vb{u) = {VK(0-; UL,u) ; ul G U}. 



Definition 2.2 A function u : [0, T] x [a, b] i—>- U is an entropy weak solution of the initial- 
boundary value problem (Ll)-(L4) on fix [0, T] x [a, 5], if it is continuous as a function from 
]0, T] into L^, and the following properties hold: 



(i) u is a distributional solution to the Cauchy problem (l.l)- f^l.^ ) on VIt in the 
sense that, for every test function 4> £ C]. with compact support contained in the set 
{{t,x) £ R^; a < X < b, t < T}, there holds 



T nb 



{u{t,x) ■ (f)t{t,x) + f{u{t,x)) ■ (j)x{t,xy)dx dt -\- I u{x) ■ (f){0,x)dx = 



(a) the flux f{u) admits weak* traces at the boundaries x — a, x = b, i.e. there exist two 
measurable functions VPa, ^'fc : [0, T] i— > R" such that 



-♦6- 



tn L°°([0, T]), 



(2.5) 
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and the boundary conditions ( l.S)-p^ are satisfied in the following sense 



*a(t) e/(Va(2aW)), t{t) & f {Vb{ub{t))) for a.e. 0<t<T; (2.6) 

(Hi) u satisfies the following entropy conditions on the decay of positive waves in time and in 
space. There exists some constant C > 0, depending only on the system so that 

(a) For any < t < T, and for a.e. a < x < y < b, there holds 



w,{t,y)^w,{t,x) <C-{^^+ log f^LJl]^ te{l,...,p}, 



(2.7) 



w,it,y)-w,it,x) <C-i^^^+\og(^^L_^jj z e {p+l,...,n}; 

(2.8) 

(b) For a.e. a < x < b, and for a.e. < ti < T2 < T, there holds 

Wi{T2,x) ~Wi{ti,x) < C - | ^^ _^^ +log (^^^1 « G {1, . . . ,_p}, 

(2.9) 

Wi{T2,x) -~w,{ti,x) < C - | ^^_^^ +log (^—^1 ^/ i e {p+l,...,n}. 

(2.10) 



Remark 2.1 The set of admissible flux values at the boundaries x = a, x = b, can be expressed 
in Riemann coordinates as 



/(Va(w)) = |/(w) ; Wi{u) = Wi{u) V i = p+ 1, . . . ,n|, 
/(Vf,(u)) = |/(w) ; Wi{u) = Wi{u) V i = 1, . . . 



(2.11) 



Hence, by the invertibility of the map f : U ^ f{U), the above boundary conditions ( ^.6| ) are 
equivalent to the set of equalities 



w^{f-\^a{t)))^w^{ua{t)) fora.c. < t < T, z = p + 1, . . . , n, 
w^{f-\^b{t))) ^w^{ub{t)) for a.e. < i < T, i = l,...,p. 



(2.12) 



This means that the boundary conditions (2.6) guarantee that, at almost every 
time t e [0, T], the solution to the Riemann problem for (1.1), having left and right initial 
states ~ Ua{t), = /^^(5'a(i)), contains only waves with negative speeds, while the solution 
to the Riemann problem with initial states = /^^('l'b(i)), — Ub{t), contains only waves 
with positive speeds. Thus, in particular, such solutions do not contain any front entering the 
domain [t, +oo[x]a,6[. 
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In the present paper we regard the boundary data as admissible controls and, in connection 
with a fixed convex, compact set T C U having the form 



r = |u e [/; Wi{u) £ [ofi, /3j], 



z = l,. 



(2.13) 



we study the basic properties of the attainable set for (L1)-(L2), i.e. of the set 
A{T) = {u{T, •) ; u 



is a sol. to (1.1) -0, Ua, U6 e L°°([0,r],r)} (2.14) 
which consists of all profiles that can be attained at a fixed time T > 0, by entropy weak solutions 



of (L1)-(L4) (according with Definition 2^) with a fixed initial data u 6 L°°([a, 6], F), and 
boundary data Ua, Ub that vary in 



(2.15) 



We will establish a characterization of ( 2.14 ) in terms of certain Oleinik type estimates on the 
decay of positive waves, and we will prove the compactness of ( p. 14 ) in the topology. 



2.2 Statements of the main results 

For any p > 0, consider the set of maps 



y- X 



KP={ 



if e L°°([a,fe], F) 



< 



y-x 



h-y 



for a.e. a < x < y < b, 
if i e {p+ I, n} 

for a.e. a < x < y < b, 
if ie{l,...,p} 



(2.16) 

The inequalities in (2.16) reflect the fact that positive waves entering through the boundaries 
x = a, X = b decay in time. Therefore, their density (expressed in terms of Riemann coordinates) 
is inversely proportional to their distance from their entrance point on the boundary. 



Theorem 2.1 Let (1.1) be a system of Temple class with all characteristic fields genuinely non- 
linear, and assume that the strict hyperbolicity condition (SHI) is verified. Then, for every fixed 
r > 0, there exists p = pir) > such that 



A{t)(ZKP yT>T. 
Moreover, there exist T > and p' < p{T), such that 

RP' C A{t) \/ t>T. 



(2.17) 



(2.18) 



Remark 2.2 Observe that, given ip G K^, any map x !—> Wi{ip{x)), i G {1, . . . , n}, is essentially 
bounded and has finite total increasing variation on subsets of [a, b] bounded away from the end 
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points a, b. Hence, any map x ^ ^^(1^9(2;)), i G {1, . . . ,71}, has also finite total variation on such 
sets and, in particular, it admits left and right limits in any point x G]a, b[. Moreover, since an 
element ip of is defined up to L""^ equivalence, we may always assume that there is a right 
continuous representative of Wi{ip), z € {1, . . . ,n}, that satisfies the inequalities appearing in the 
definition of K^. 



Theorem 2.2 Under the same assumptions of Theorem 1, the set A(T) is a compact subset 
ofL^{[a,b], r) for each T > 0. 

Indeed, we will prove in Section ^ that the compactness of the attainable set A{T) holds even 
in the case where A{T) is defined as the set of all configurations that can be reached at time T 
only by solutions of the mixed problem for (1.1) that admit a strong trace at the boundaries 
X — a, X = b (as the ones generated by a front tracking algorithm). 



2.3 Existence and uniqueness of solutions 



We describe here a front tracking algorithm that generates approximate solutions to (1.1) on 
the strip [a, b] x R"*" continuously depending on the initial and boundary data, which represents 
a natural extension of H, |l2j . Fix an integer v > 1 and consider the discrete set of points in F 
whose coordinates are integer multiples of 2^"^: 

F'' = |u e F; Wi(u) e 2-''Z, i = l,...,7i|. (2.19) 

Moreover, consider the domain 

2?''^|(u,m',u"); u e L°°([a,6],F''), -u',u" G L°°(R+,F''), u, u' , u" are piecewise constant |. 

(2.20) 

On we now construct a flow map whose trajectories are front tracking approximate solu- 
tions of (Fl). To this end, we first describe how to solve a Riemann problem with left and right 
initial states u^,u^ G F''. In Riemann coordinates, assume that 

f L\ ■ L r L L\ r R\ ■ R. I R R\ 

Consider the intermediate states 



z 



z'-ii(i«f,...,i.f, u;^), z"=u«. (2.21) 



The solution to the Riemann problem (u^, u^) is constructed by piecing together the solutions to 
the simple Riemann problems (z*~^,z*), i = 1, . . . ,n. If wf^ < wf', the solution of the Riemann 
problems {z'^~^, z^) will contain a single i-shock, connecting the states z'"^, z*, and traveling with 
the Rankine-Hugoniot speed Ai(z*~^, z'). Here and in the sequel, by Xi{u,u') we denote the i-th 
eigenvalue of the averaged matrix 

A{u,u')=[ Df{eu + {i-e)u')de. (2.22) 
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If wf^ > , the exact solution of the Ricmann problem (z^~^,z^) would contain a centered 
rarefaction wave. This is approximated by a rarefaction fan as follows. If wf = wf + pi 2^'^ we 
insert the states 

z'''^{w^,...,wf+i2-'',wf+,,...,w^), e = 0,...,p,, (2.23) 

so that z*'" = z^'P^ = z*. Our front tracking solution will then contain pi fronts of the i-th. 

family, each connecting a couple of states z*'^ and traveling with speed Ai(z*'^~^, z*'^). 

For any given triple of (piecewise constant) initial and boundary data {u, Ua , Uf,) € V^, the 
approximate solution u{t, •) = E^{u, Ua , Ub) is now constructed as follows. At time t = 0, for 
a < X < b we solve the initial Riemann problems determined by the jumps in u according to the 
above procedure, while at a; = a we construct the solution to the Riemann problem with left and 
right initial states — Ua{0+), — u{a+) and take its restriction to the interior of the domain 
Q. In the same way, at x — b, we take the restriction to the interior of il of the solution to the 
Riemann problem with initial states = u{b—), = Ub(0+). This yields a piecewise constant 
function with finitely many fronts, traveling with constant speeds. The solution is then prolonged 
up to the first time ti at which one of the following events takes place: 

a) two or more discontinuities interact in the interior of Q; 

b) one or more discontinuities hit the boundary of Q] 

c) the boundary data Ua has a jump; 

d) the boundary data Ub has a jump. 

If the case a) occurs, we then solve the resulting Riemann problems applying again the above 
procedure, while in the other three cases b)-c)-d) we construct the solution to the Riemann 
problem with left and right initial states — Ua{ti+), = u{ti,a+), or — u{ti,b~), 
= 'Ub{ti+), and take its restriction to the interior of the domain fl. This determines the 
solution u{t, •) until the time t2 > ti where one of the events a),b),c) again takes place, etc. . . 
Notice that at any time where case b) occurs but c) or d) do not take place, no new wave is 
generated. Therefore, waves entering the domain Q at the boundaries x ~ a, x ^ b are produced 
only by the jumps of the boundary data Ua, Ub- 

As in 1^, ^] , one checks that the approximate solution u constructed with this algorithm is 
well defined for all times t > 0. Indeed, the following properties hold. 

- The total variation of u{t, •), measured w.r.t. the Riemann coordinates wi{t, ■),... , Wn{t, •), 
is non-increasing in time. 

- The number of wave-fronts in u(t, ■) is non-increasing at each interaction. Hence, the total 
number of wave-fronts in u{t, •) remains finite. 

It is then possible to define a fiow map 

p ^ E1y>. p = {u,Ua,Ub)€V, t>0 (2.24) 
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of approximate solutions of (pTl|). By construction, each trajectory t i— > E^p is a weak solution 
of (1.1) (because all fronts of u{t,-) = E^p satisfy the Rankine-Hugoniot conditions), but may 
contain discontinuities that do not satisfy the usual Lax stability conditions (due to the presence 
of rarefaction fronts). On the other hand, one can verify as in Lemma 4.4] that, due to genuine 
nonlinearity, the amount of positive waves in u{t, •), measured w.r.t. the Riemann coordinates 
wi(t, •),... , Wn{t, •), decays in time and in space. Hence, for a.e. a < x < y < b, one obtains the 
Oleinik type estimates 



w,{t,y)-w,{t,x)< C-iy—^+\og(^LJl\K+N,2-'' if ie{l,...,p}, 



t \x — b 

w^{t,y)-w^{t,x)< C-{y-^+\og(^^^]\+N,2-'' if z G {p + 1, . . . , n}. 



(2.25) 



where Ni, denotes the maximum number of shocks of each family present in the initial data u, 
and in the boundary data Ua , Ub- Similarly, one can check that along the x-sections, for a.e. 
< Ti < r2, there holds 



W^{T2,x)-W^{Tl,x)< C-{^-^^^^—^+\0g(^^^]\+N,2-'' if ie{l,...,p}, 

w^{t2,x)-w.,{ti,x)< C-{^^—^+\og(-]\+N,2-'' if i e {p+l,...,n}. 

X ^ a \Ti 



(2.26) 



Remark 2.3 Observe that, if u{t,x) is a front tracking solution of the Cauchy problem for (1.1) 



(with initial data u{x) = u{0,x)) constructed by the algorithm in |12 on the upper half 
plane M+ x M, then the restriction of u{t, •) to the interval [a, b] coincides with the front tracking 
solution E^{u, Ua,Ub) of the mixed problem for (1.1), with boundary data Ua{t) u{t,a), 
Ub(t) = u{t, b) . 



As — > oo, the domains "V^ become dense in 

V=^{u,UaUb) ; u e L°°([a,5],r), Ua, £ L°°(R+, T)}. 



(2.27) 



Thus, following the same technique adopted in B, one can define a flow map Et on P as a suitable 



limit of the flows E^ in ( 2.24 ), that depends Lispschitz continuously on the initial and boundary 
data. Namely, the following holds. 



Theorem 2.3 Let (1.1) be a system of Temple class with all characteristic fields genuinely non- 
linear, and assume that the strict hyperbolicity condition (SHI) holds. Then, there exists a 
continuous map 



{t, U, Ua, Ub) ^ Et{ U, Ua,Ub) 



t>o, (u, Ua, Ub) ev , 



(2.28) 



and some constant C > depending only on the system (1-1) and on the domain T, so that. 
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for every fixed < S < {b — a)/2, and for all pi = (u, Mq , Ub ), P2 = (v, Va , Vb) € T>, letting 
Lt = Lt{5) ^ C(l + \og{t/5)), there holds 

WEtPl - ^tP2||Li([a+5.6_5]) ^ ■ " - + II " II ([0, *]) + II " II ([0, t]) } 

(2.29) 

for all t > S. Moreover, the map {t,x) <—>■ Et{u, Ua , Ub){x) yields an entropy weak solution (in 
the sense of Definition 2.i) to the initial-boundary value problem ( Ll)- ^!^ ) on Q,, that admits 
strong traces at the boundaries x — a and x ~ b, i.e. there exist two measurable maps 
ipa, ipb ■ ^ U .such that 



lim / \Et{u, Ua , Ub)ix) - ^a{t)\ dt = 0, 



lim 



V r > 0. 



(2.30) 



\Et(^U, Ua , Ub){x) - 'ljjbit)\ dt = 0, 



The proof of Theorem p.3| can be obtained with entirely similar arguments to those used to 



establish Theorem 2.1], where a continuous flow of solutions to (1.1) is constructed in the case 
of a mixed problem on the quarter of plane {(t, x) G ; t > 0, x > 0}, with a single boundary 
at x = 0. 



Concerning uniqueness, with the same arguments in Q one obtains the following result which 
is the extension of Theorem 2.2] to the present case of a domain Q with two boundaries at 
X = a and at a; = 6. 



Theorem 2.4 Let (1.1) be a system of Temple class satisfying the same assumptions as in The- 



orem 2.S. Let u = u{t,x) be an entropy weak solution to the mixed problem (l.l)-(1.4) on the 
region Ht == [0, T] x [a,b] (in the sense of Definition \2.^ . Assume that the following conditions 
hold. 



(i) The map {t,x) (u{t,-), u{-,x)) takes values within the domain 

Vt ^ {{u, Ua, Ub) ; ueL°°([a,6],r),S„«6eL°°([0,T],r)}. 

(a) There holds 



ess sup / \u{t,x) — u{x)\dx — Q . 



(Hi) There holds 



ess sup / \wi(^u(t, x)) — Wi{ua{t))\ dt — V i = p + 1, . . . , n, 

x^a+ Jo 



ess sup 



(u{t, a;)) — Wi{ub{t))\ dt — \/ i — 1, . . . ,p. 



(2.31) 
(2.32) 

(2.33) 
(2.34) 
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Then, u coincides with the corresponding trajectory of the flow map Et provided by Theorem \2.^ , 
namely one has 

u{t,-)^ Et{u,ua,ui,){-), yO<t<T. (2.35) 



The next result shows that the conditions ( ^.32| )-( ^r34 ) are certainly satisfied by entropy weak 
solutions to the mixed problem ( |l.l| )-( pr^ ) obtained as limit of front tracking approximations. 



Theorem 2.5 Let (1.1) be a system of Temple class satisfying the same assumptions as 
in Theorem 2.S. Consider a sequence ■) : [a,b] F" of wave-front tracking approximate 

solutions of the mixed problem for ( LA ) ( constructed with the above algorithm ) that converges 
in Ti^ , as V oo, to some function u{t,-) : [a,b] i-^ F, for every t G [0,T]. Then, there exist 
the right limit at x — a, and the left limit at x = b, of the map x — ^ u{t,x) for every t £ [0,T], 
and the right limit at t = of the map t u{t,x) for every x G [a.,b]. Moreover, there is a 
countable set A/" C M such that u{t,a) — u(t,a+), u(t,b) — u{t,b—) for all t G [0,T] \ TV, and 
u{0, x) = u{0+, x) for all x G [a, b] \ M , and, setting u = u{0, ■), Ua a), = h), there 

holds ^Idj ). 



Remark 2.4 It was shown in P, Lemma 2.1] that an alternative way to prove the essential 



limits (2.33)-(2.34), is to employ the distributional entropy inequalities associated to the "bound- 
ary entropy pairs" for (p^), introduced by G.-Q. Chen and H. Frid in ||lj, |l5|. However, in order 
to apply 1^, Lemma 2.1] to a function u obtained as a limit of approximate solutions u'^, it is 
necessary to know the convergence of the sequence of the corresponding boundary data u'^, u'^. 
Instead, the result provided here by Theorem 2.5 allows to derive the limits ( ^.33 )-(2.34) requir- 
ing only the convergence of the sequence of the approximate solutions u'^{t, ■), for all t. This 
property will be crucial to establish the main result of the paper stated in Theorems p.l|-^.2|. 



In order to prove Theorem 2.5, we will show in the next section that, for Temple systems, 



solutions of the mixed problem (1.1)- (1.4) with possibly unbounded variation enjoy the same 
regularity property (of being continuous outside a countable number of Lipschitz curves) possessed 
by solutions with small total variation of a general system, thus extending the regularity results 
obtained under the smallness assumption of the total variation by DiPerna [|l^ and Liu 
(for solutions constructed by the Glimm scheme) and by Bressan and LeFloch ]|l^ (for solutions 
generated by a front tracking algorithm). 



Proposition 2.1 In the same setting as Theorem 2.5, consider a sequence u'^(t,-) : [a,b] i— * F"^ 
of wave-front tracking approximate solutions of the mixed problem for j| j ( constructed with the 
above algorithm) that converges in , as v ^ oo, to some function uit,-) : [a,b] F, for every 
t G [0,T]. Then, there exist a countable set of interaction points Q = {{ti, xi); I G N} C J^t = 
[0, T] X [a, b], and a countable family of Lipschitz continuous shock curves T = {a; = ym{t); t G 
]^rm Sjyii, fn G N}, such that the following hold. 
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(i) For each ni G N, and for any r S]r„i, s„j[ with {T,ym{T)) ^ Q, there exist the derivative 
ym{T) and the left and right limits 

lim u{s,y)^u-, lim u{s,y) = u+ . (2.36) 

{s,y)~*{T,y„^{T)),y<y^{T) (s,a)^(T,a„, (t)) , j/>y„ (r) 

Moreover, these limits satisfy the Rankine Hugoniot relations 

ym{T)-{u+ -u-) = f{u+)- f{u-) (2.37) 
and, for some i G {1, ...,n}, there hold the Lax entropy inequalities 

K{u+) < y,n{t) < X^{u-) . (2.38) 

(a) The map u is continuous outside the set U T. 



3 Preliminary results 

In this section we first provide some estimates on the distance between two rarefaction fronts 



of a front tracking solution (constructed by the algorithm described in Section 2.3) similar 
to [|2[ Lemma 4], |0, Prop. 4.5]. We next show how to approximate the profile u{t,-) of a 
solution of the mixed problem (pTl|)-(pr^), with a function taking values in the discrete set 



defined at (|2.19| ), which enjoys the same type of estimates on the positive waves as u{t, ■). We 
conclude the section establishing the regularity result stated in Proposition on the global 
structure of solutions to the mixed problem for (^j^), which in turn yields Theorem p^ . 



Lemma 3.1 There exists some constant Ci > depending only on the system (1.1) such that 
the following holds. Consider a front tracking solution u{t,x) with values in , constructed by 
the algorithm of Section \2.^ on the region [r, r'] x [a,h]. Then, given any two adjacent rarefaction 



fronts of u located at x{t) < y{t), t G [r, t'], and belonging to the same family, there holds 

\y{r') - xiT')\ < \y{T)-xiT)\+CiiT' -t)2-- . (3.1) 

Proof. Consider two adjacent rarefaction fronts of the fc-th family x{t) < y{t), t G [r, r'], and 
let Ti < ... < TN be the interaction times of x{t) in the interval [r, r']. Set tq = r, tn+i = t' , 
and fix a G {0, N}. Let t z{t; s, x) be the characteristic curve of the fc-th family starting at 
{s,x), i.e. the solution to the ODE 

z = Xk{u{t, z)), z{s;s,x)=x. 

Notice that, although the above ODE has discontinuous right hand-side (because of the discon- 
tinuities in the front tacking solution u), its solution z{-;s,x) is unique and depends Lipschitz 
continuously on the initial data x since it crosses only a finite number of jumps (see 0). Choose 
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to < ti < Ta+i so that the characteristic curve z{-;to, x{t(j)) does not cross any wave-front of 
the other famihes in the interval [to, ti], and then, by induction, define a sequence of times 
{tijiez C ]tq,Tq+i[ so that 



Ta < t-i-l < t-t <to<ti< ti+i < Ta + l, i G N, 



hni ti =Ta, hm ti = r^+i, 

— » — oo i — >+oo 



(3.2) 



with the properties that the characteristic curve of the fc-th family starting at (ti, x{tiy), does not 
cross any wave-front of the other families in the interval \ti, ti+i], for each i ^"L. Thus, setting 

ul = u{ti, x{ti)+), u'l = u{ti, x{ti)-), 

and observing that, by construction, one has \w{u'^) — w{u^)\ < 2^'^ , we derive 

\z{ti+i; ti, x{ti)) - x{ti+i)\ < (ti+i - ti) ■ |Afe(M+) - Afc(u+, u^)\ 



< c ■ {ti+i - ti) ■ |w(u+) - )| 



(3.3) 



for some constant c > depending only on the system. Relying on (3.3), and since z(t'; ii+i, x) 
depends Lipschits continuously on the initial data x, we deduce that there exists some other 
constant c' > 0, depending only on the system and on the set F, so that there holds 



\z{t'; ti, x{ti)) - z{t'; ti+i, x{ti+i))\ < c' ■ \z{ti+i; ti, x{ti)) - x{ti+i)\ 

<c' -c- {t,+i - t{) ■ 2-" 



(3.4) 



for any i £ Z. Thus, by (3.2), and thanks to ( |3.4D , we obtain 

|z(t'; Ta, x{Ta)) - z{t'; T^+i, x{Ta))\ < ^ \z{t'; ti,x{ti)) - z{t' ; ti+i, ,x{ti+i))\ 



(3.5) 



< c' • c • (r„+i - T„) • 2-'^ . 
Repeating this computation for every interval ]tq,Tq-|_i[, a € {0, ...,N}, we get 



N 



\z{t'; T, x{t)) - x{t')\ < ^ \z{t'; Ta, x{Ta)) - z{t'; Ta + l, x{Ta))\ 



<c'-c-{t'-t)- 2-"" . 



(3.6) 



Clearly, one obtains the same type of estimate as ( pig ) for the other rarefaction front y{t), i.e. 
there holds 

\z{t'; t, y(r)) - y(r')| < c' ■ c ■ {t' - t) ■ 2-" . (3.7) 



On the other hand, by (2^), we have 

|z(t'; r, x(t)) - z{t'; t, y(r))| < \x{t) - y(r)| -|- 2 A"-'^ • (r' - r). 



Thus, rt3 



together yield (3.1), concluding the proof. 



(3.8) 
□ 
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In the following, in connection with any (right continuous) piecewise constant map 
Tp : [a, b] 2^'^ Z, we will let tt^tJj) — {xq — a < xi < ■ ■ ■ < xj = b} denote the partition 
of [a,b] induced by ip, in the sense that 'ipi^) is constant on every interval [xg, < £ <£. 

Then, given p > 0, for any v > 1, consider the set of piecewise constant maps 



Wi{(p{xk)) - Wi{ip{xh)) 



< 



bp 



Xh - a 



ip:[a,b]^V^; 



Wi{ip{xk)) - Wi{ip{xh)) 



< 



5p 



Xk ~ Xh 



b-Xk 



for 



if 



for 



if 



a < Xh < Xk <b, 
Xh-, Xk e 7r(wi o ip), 

i e {p+ 1, ...,n} 

a < Xh < Xk < b, 

Xh, Xk e TT{Wi O ip), 

i e {l,...,p} 



(3.9) 

The next lemma shows that we can approximate in any map (p G with a piecewise constant 
function ip^ £ K^. 



Lemma 3.2 For any given p G , there exists a sequence of right continuous maps ip^, £ K^, 
V >!, such that: 

a) for every i £ {1, . . . , n}, and for any Xh £ T^{wi o Lpu), there holds 

Wi{p}v{xh+i)) > Wi{pu{xh)) Wi{p{xh+i)) ^ Wi{p{xh)) + 2"'' ; (3.10) 



b) there holds 

p^^^p in Li([a,6]). (3.11) 

1. First observe that, by Remark ^.2| , any map x ^ Wiikpixf), i £ {1, . . . , n} has finite total varia- 
tion on [a + e, 6 — e], £ > 0. Hence, we may assume that Wi(ip(/)) admits left and right limits in any 
point x £]a, 6[, and that Wi(}p{x)) = Wi{(p{x^)) = lim^^j,+ Wi{p{^)), for all i £ {1, . . . ,n}. Let 
{yijn ; fTT- G N} be the countable set of discontinuities of Wi{ip{-)), i £ {l,...,n}. Then, we 
can find a partition ^^^^ ^ yi^^n < < ■ ■ < = Vi.™-' of each interval yi^m\ where 

x I— > Wi{}p{x)) is continuous, so that: 

i) for every \ < £ < £i,„i there holds 

w.{p{^lJ)e2-^Z; (3.12) 



a) for every 1 < £ < £i,m one has 

\w,{p^{x)) - w,{p,{il„:))\ <2-^ V a; £ ^^[. (3.13) 
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Notice that the Oleinik type conditions stated in the definition of if imply that, at any discon- 
tinuity point yi^m of Wi{ip{-)), one has 

hm w,{(f {£,))> w,{ip{yi^,n)) ■ (3.14) 



2. Let ip^ : [a, b] i— s- be the piecewise constant, right continuous map defined by setting, for 
every i e {1, . . . , n}, and for any interval [yi^rm l/i.m' [ where Wi{ip{-)) is continuous. 



2-l2-u;.(^(ei™))J if 
2--(L2-z«,(^(ei„))J+l) if 



X ^[^Ln^ &,m[^ and 
^^[^Ln^^lmi^ and 

w^{M^lJ) > 2--(L2-u;.Mei„J)J + 2-i), 



(3.15) 

where [-J denotes the integer part. Notice that, by construction, and because of (3.12)-(3.13), 
(3.14), the map ipi, : [a, b] i— > enjoys the following property 

Xh < Xk € TT{Wi O ip^) 



Wi{(p{xk)) > Wi(ip{xh)) + 2 



(3.16) 



Therefore, since G K^, relying on ( 3.13| ), ( 3.16| ), we deduce that, for every Wi{Lpi,{-)), 
i e {p + 1, . . . ,n}, and for any Xh < Xk & TT{wi o ip^) such that Wiiip^ixk)) > Wi{ipi,{xh)), 
there holds 



Wi(ipt,ixk)) - Wi{ipi,{xh)) ^ Wi{ip{xk)) - Wt{ip{xh)) + 2 ^) 



< 



5 {wi{ip{xk)) - Wi{ip{xh))) 



Xk - Xh 



< 



5p 



Xh - a 

Clearly, with the same computations, we can show that, for every Wi((p^(-)), i G {1, 
for any Xh < Xk ^ Tr{wi o (pu), there holds 



Wi{ipi,{xk)) - Wi{ipy(xh)) 



< 



Xk - Xh 



5p_ 

b- Xk 



(3.17) 



and 



(3.18) 



The estimates (3.17)-(3.18), together, imply that G K^, while (3.13) yields (3.11). On the 
other hand observe that, by construction, and because of (3.14), the map (p,y satisfies condition 
(3.10), which completes the proof of the lemma. □ 



We now provide a further estimate on the distance between two rarefaction fronts of a front 
tracking solution that, at a fixed time r, attains a profile belonging to the set (3.9). 
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Lemma 3.3 Consider a front tracking solution u{t,x) with values in , v > 1 , constructed by 



the algorithm of Section 2.S on the region [r, t'] x [a, 6]. Assume that u(t' , ■) is right- continuous, 



verifies condition a) of Lemma 5.1, and satisfies 

uiT',-)GKP\ 



P 



6Ci 



(3.19) 



where A™™, Ci, are the minimum speed in (2.5), and the constant of Lemma 3.1. Then, given 
any two adjacent rarefaction fronts of u located at x(t) < y{t), t G [r, r'], and belonging to the 
same family, there holds 

x{t) < v{t). (3.20) 
Proof. To fix the ideas, assume that x{t) < y{t) are the locations of two adjacent rarefaction 



fronts of the fc G {p + 1, . . . , n} - th family, and hence, by (2.2), have positive speeds. Observe 



that, by condition a) of Lemma 3.2, one has 



Wk{u{T', y{T'))) - wu{u{t', x{t'))) = 2-'' . 



(3.21) 



Moreover, since u is a front tracking solution constructed by the algorithm of Section 2^ on the 



region [t, t'] x [a, 6], we can apply Lemma Thus, using ( p.3|) , (3J), ( 3.21 ), and recalling the 



definition (3J) of , we deduce 



y(r') - x{t') < y{T) - x{t) + Ci(r' - r) 2"^ 

x{t') — x{t) 



<y{r)-x{T)+Ci- 
< y{T) ~ x(t) + Ci 



A" 



5p' 



(y(r')-x(r')) 



which, because of (3.19), implies 

yir) - x{t) >{l- • W) - > . 



proving (3.20) 



□ 



We next derive a regularity property enjoyed by general BV solutions of Temple systems 



defined as limit of front tracking approximations, which allows us to establish Proposition 2.1 



This is an extension of the regularity results obtained in [0 ^ ^ for solution with small total 
variation of general systems. The arguments of the proof are quite similar as for the corresponding 
result in ||l^, but we will repeat some of them for completeness, referring to ||l^ (see also [|[ 
Theorem 10.4]) for further details. 

Lemma 3.4 Let ^l.f\ ) be a system of Temple class satisfying the same assumptions as 



in Theorem 2.c. Consider a seguence u'^{t,-) : [c, d] ^ P", t G [r, s], of front tracking ap- 



proximate solutions of the mixed problem for (1.1) (constructed by the algorithm of Section 2.c), 
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that converges inlJ- , as v oo, to some function u{t, •) : [c, d] T, for every t € [r, s] C M"*". 
Assume that 

Tot.Var.{u''{t,-)) < M, Tot.Var.iu" {■ ,x)) <M M t,x,v, (3.22) 

for some constant M > 0. Then, there exist a countable set of interaction points 
Q = {(t;, xi)] I G N} C -D = [r, s] X [c, d\, and a countable family of Lipschitz continuous 
shock curves T = {x = ym{t)', t £ ]^m, Sm[, rn dz N}, such that the following hold. 

(i) For each m G N, and for any r G]r„i, Sm[ with {T,ym{T)) ^ Q, there exist the left and 
right limits ( 2.3t ) of u at {T,ym{T)) and the shock speed ym{T). Moreover, these limits 
satisfy the Rankine Hugoniot relations (2.37) and the Lax entropy inequality (2.38), for 
some i G {1, n}. 

(ii) The map u is continuous outside the set U T. 

Proof. 

1. To establish (i) we need to recall some technical tools introduced in (see also The- 
orem 10.4]). For every front tracking solution u"^, we define the interaction and cancellation 
measure fil^ that is a positive, purely atomic measure on D, concentrated on the set of points P 
where two or more wave- fronts of u'^ interact. Namely, if the incoming fronts at P have size 
ai, . . . ,ae (w.r.t. the Riemann coordinates), and belong to the families ii, . . . , respectively, we 
set 



{ia ; ia — 



(3.23) 



Since ^If^ have a uniformly bounded total mass, by possibly taking a subsequence we can assume 
the weak convergence 

/if - (3.24) 

for some positive, purely atomic measure /Lt^*^ on D. Call Q the countable set of atoms of fJ,^'~^ , 
i.e. set 

e = {P G D; /^(F) > 0}. 

For every approximate solution u'^ taking values in F", > 1, and for any fixed e > 2~'^, 
by an e~ shock front of the i—th family in u'^ we mean a polygonal line in D, with nodes 
(tq, xo), (tjv, xn), having the following properties. 

(I) The nodes {Th,Xh) are interaction points or lie on the boundary of D, and the sequence of 
times is increasing tq < ti < ■ ■ ■ < tn . 

(II) Along each segment joining (r^^i, Xh-i) with (r^, Xh), the function u'^ has an i— shock with 
strength \ah\ > e. 

(Ill) For h < if two (or more) incoming i— shocks of strength > e interact at the node (rh, x^), 
then the shock coming from {rh-i, x^^i) has the larger speed, i.e. is the one coming from 
the left. 
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An e— shock front which is maximal with respect to the set theoretical inclusion will be called 
a maximal e— shock front. Observe that, because of (III), two maximal e— shock fronts of the 



same family either are disjoint or coincide. Moreover, by (3.22), the number of maximal e— shock 
front that starts at the boundary of D is uniformly bounded by 3M/e. On the other hand, the 
special geometric features of Temple class systems guarantee that no new shock front can arise 
in the interior of D. Indeed, the coinciding shock and rarefaction assumption together with the 
existence of Riemann invariants prevents the creation of shocks of other families than the ones of 
the incoming fronts at any interaction point. Therefore, for fixed e > 0, and i G {1, . . . ,n}, the 
number of maximal e— shock front of the i-th family remains uniformly bounded by AI^ = 3M/e 
in all u'^, J/ > 1. Denote such curves by 

ylm ■■ K;™, it;™] ^K, m = i,...,m,. 

By possibly extracting a further subsequence, we can assume the convergence 
J/t,m(-) ^2/™(-), tl%^t'^^, m^l,...,M,, 

for some Lipschitz continuous paths : [i^", tf^'^] i-^ R, m = 1, AI^. Repeating this construc- 
tion in connection with a sequence et — > 0, and taking the union of all the paths thus obtained, 
we find, for each characteristic family i G {1, ... ,71} , a countable family of Lipschitz continuous 
curves j/m : [t^, t+] i-^R, m G N. Call T the union of all such curves. 

2. Consider now a point P = (r, ym^r)) ^ O along a curve y„i G T of a family i G {1, . . . , n}. 



Notice that, by construction, and because of (3.24), no curve in T can cross at P. Moreover 



by ( |3.22 ), the function u{t, •) has bounded variation, and hence there exist the limits 



lim u{t,x)=:u , lim u{t,x)^u'^. (3.25) 



We claim that also the limits (2.36) exist, and thus coincide with those in ( [3.25| ). To this end 
observe that, by construction, there exist a sequence of shocks curves yu,m of the i-th family 
converging to y™, along which each approximate solution u'^ has a jump of strength > £*, for 
some e* > 0. Then, relying on the assumption 

/^({P}) = 0, (3.26) 

and letting B{P,r) denote the ball centerd at P with radius r, one can establish the limits 

lim limsup sup lu"^ (t, x) — =0, (3.27) 

""^0+ 1^^+00 y (t,x)e-B(P,r), x<y„,„(t) J 

lim hmsup sup \u'^{t,x) — =0, (3.28) 



which clearly yield ( ^.36| ). Indeed, if for example ( 3.27 ) do not hold, by possibly taking a subse 



quence we would find e > and points Pi, = {t,y, — > P on the left of y^,m such that 
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On the other hand, by the first hmit in ( |3.25 ), and since u'^{t,x) — > u{t,x) for a.e. x £ [o:,f3], we 
could also find points = (r, ^ P on the left of y^,m such that 



> A" 



where X^^'^ denotes the maximum speed at (p.3[). But then, for each solution u'^ , the segment 



-Py Qv would be crossed by an amount of waves of strength > s. Hence, by strict hyperbolicity 
and genuine nonlinearity, this would generate a uniformly positive amount of interaction and 
cancellation within an arbitrary small neighborhood of P (see. [|[ Theorem 10.4-Step 5]) which, 
by the definition ( 3.23| ), and because of ( ^.24| ), contradicts the assumption ( 3.26| ). 



To complete the proof of (i) observe that, by construction, the states u~^(t), w+^(t) to the 
left and to the right of the jump in at yu,m{T) satisfy the Rankine Hugoniot conditions. Thus, 
relying on ( ^.27 )-( 3.2g ), and on the convergence yi^^m Vu, one deduces ( 2.371 ). The proof of (ii) 
can be established with the same type of arguments (cfr. ||, Theorem 10.4-Step 8]). □ 



As an immediate consequence of Lemma 3.4, we derive Proposition 2.1 stated in Section 2.3 



Proof of Proposition |2.l| . Consider a sequence u''{t, ■) : [a, b] ^ of front tracking approx- 



imate solutions of the mixed problem for (1.1) on the region JIt == [0, P] x [a, 6], that converges 
in L^, as — > cxd, to some function uit, ■) : [a, b] > P, for every t e [0, T]. Observe that, by The- 
orem 2.3 one can find another sequence {u'^},y>i of approximate solutions of (1.1) on the region 



fi-r, whose initial and boundary data have a number of shocks < v for each characteristic 
family, and such that 



l"'^(^.-)-«'^(i,-)||Li([„,,])<l/^ 



V t e [i/u, T] . 



Then, thanks to the Oleinik estimates (2.25)-(2.26), and because all take values in the compact 



set (2.13), there will be, for every fixed e > 0, some constant > such that 

Tot.Var.{u''(t, ■) \ [a + e, b - e]} < V < e [e, T] , 

Tot. Var. {«''(■, x) ; [e, T]} < V a; G [a e, 6 - e] , 

Thus, writing as the countable union 

= i^kDk, Dk = [l/k, T]x[a+ {\/k), b - (1/fc)] 



V e N. 



(3.29) 



and applying Lemma 3.4 to each sequence of maps vY = v'^ Id^: ^ ^ li defined as the restriction 



of v'^ to the domain D^, we clearly reach the conclusion of Proposition 2.1 



□ 



We are now in the position to establish Theorem 2.5, relying on Proposition p.l| and on The 



orem 2.4 



Proof of Theorem p.5| . Let uy(t, •) : [a, 6] i-^ P'' be a sequence of front tracking approximate 
solutions of the mixed problem for (|l.lD on the region ilj- = [0, T\ x [a, 6], that converges in L^, 
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as — > oo, to some function •) : [a,b] t— > F, for every t G [0,T]. Since, by construction, 
each u"^ is a weak solution of (pTl|), and because ^"^(0, •) — * u(0, •) = u, also the limit function 
u is a weak solution of the Cauchy problem (1.1)-(1.2) on the region J^t- Moreover, applying 
Proposition 2J, we deduce that u admits at i = and at a; = a, a; = 6 the left and right 
limits stated in Theorem 2.5. On the other hand, by the same arguments used in the proof 
of Proposition 2A, we may assume that the initial and boundary data of each approximate 
solution u'^ have at most < v shocks for every characteristic family. Then, letting v ^ co 
in ( ^.25| )-(2.26), by the lower semicontinuity of the total variation we find that u satisfies the 
entropy conditions (2.7)-(2.1C) on the decay of positive waves. It follows that u is an entropy 
weak solution of the mixed problem ( |l.lD -(1.4) according with Definition 2.2. Hence, observing 
that by constructin the map (t, x) {u{t, •), x)) takes values within the domain Vt defined 
in ( ^.31 ), and applying Theorem |2.4| , we deduce that ( 2.35| ) is verified. □ 



4 Proof of Theorems g^-^ 



Proof of Theorem 2.1. We shall first prove that, for every fixed r > 0, there exists some 



constant p = p{t) > so that ( p.l7| ) holds. Given Ua £ U^, ut G U^, r > r, let u = u{t, x) be an 



entropy weak solution of (1.1)-(1.4) on the region [0, t] x [a, b] according with Definition 2.2. Then, 
the Oleinik-type estimates (|2.8| ) on the decay of positive waves imply that, for i e {p + 1, n}, 
r > r, and for a.e. a < x < y < b, there holds 



w»(t, y) - Wi{T,x) , ] 

y-x ~ 



log 



<{b-a)C ■ {- + 

T X — a 



X — a 
1 



(4.1) 



< 



C{b-a){{b-a) + T) 



X — a 



Clearly, with the same computations, relying on the Oleinik-type estimates (2.7), we deduce that, 
for i e {1, T > T, and for a.e. a < x < y < b, there holds 



Wi{T,y) - W^{T,X) ^ g(6-a)((6-Q)+T) 



y-x 



Hence, taking 



P> 



C(6-a)((5-a)+T) 



(4.2) 



(4.3) 



from (4.1)-(4.2) we derive u(r, •) £ ^ which proves (2.17). 



Concerning the second statement of the theorem, we will show that, letting A™'", p', be the 
minimum speed in ( ^.3| ), and the constant ( |3.19 ) of Lemma 3.1, and taking 



T 



4 (6 - g) 



(4.4) 
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the relation ( 2.1§| ) is verified, i.e. that, given ip G , and t > T, there exist Ua G 



Uf, e U!^, and a solution u{t,x) of (1.1)-(1.4) on [0, r] x [a, 6] (according with Definition 2.2), 
such that w(t, •) = ip. Notice that, by Remark |2.2| , we may assume that Wi((p(x)) admits left and 
right limits in any point x G]a, 6[, and that Wi{ip{x)) = Wi(ip(x^)) = lim^^^+ 'Wi{ip{£_)), for all 
i e {1, . . . , n}. The proof is devided in two steps. 

Step 1. Backward construction of front tracking approximations. Letting p' > be the 
constant in ( 3.19| ), consider a sequence {^u}u>i of (right continuous) piecewise constant maps 



in , satisfying the conditions a)-b) of Lemma 3.2, and take a piecewise constant approximation 



of the initial data u, so that 



u in L^. Given t > T (T being the time 



defined in ( |4.4D ), for each v > 1, we will construct here a front tracking solution u'^{t, x) of (1.1) 
on the region [0, t] x [a, 6], with initial data ^"^(0, •) — u'^, so that 

u^T,-)=cp,. (4.5) 

This goal is accomplished by proving the following two lemmas. 

Lemma 4.1 Let T, p' > he the constants in {4ji) and ( 3.1!\ ). Then, for every (right continu- 



ous) £ RP , V > 1, satisfying the condition a) of Lemma 3.i, and for any t > T, there exists 
a front tracking solution u'^{t,x) of ( [_?. on the region [(3/4)T, t] x [a, 6], with boundary data 
51^ ^ u^i-, a), = u^i-, b) e L°°([(3/4)r, r], F'^), so that 



w''((3/4)T, x)=L,, 
for some constant state G F". 



'(r, x) = (p^{x) , 



V a; G [a, 6] 



(4.6) 



Proof. Given t > T, and ifi^ G KP , v > 1, satisfying the condition a) of Lemma 3.2, we will 



use the algorithm described in Section 2.3 to construct backward in time a front tracking solution 
that takes value ip^ at time r. To this end, we first observe that according with the algorithm of 



Section 2.3 we can always construct the backward solution of a Riemann problem with terminal 
data 



u{t, x) = 



if 
if 



(4.7) 



if the the terminal states u^, G F"^ have Riemann coordinates 



w 



that satisfy 



Indeed, if we consider the intermediate states 

,L 



u 
,,R 



■ I ^n—t I 



W. 



R 

n— i+1 ' 



if 
if 
if 



Vi. 



i = 0, 

< Z < 71 , 

1 = n , 



(4.8) 



(4.9) 
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we realize that, because of (4.8), the solution of every Riemann problem with initial 



states (z* , z*) (defined as in Section |2.3D contains only a single front. Thus, we can construct the 



solution to the Riemann problem with terminal data (4.7) in a backward 
neighborhood of (i, ^) by piecing together the solutions to the simple Riemann problems (z*~^, z'), 
i = 1, . . . , n. 

A front tracking solution can now be constructed backward in time starting at i = r, and 
piecing together the backward solutions of the Riemann problems determined by the jumps in ip,y. 
The resulting piecewise constant function u'^{t—, •) is then prolonged for t < t tracing backward 
the incoming fronts at t = r, up to the first time ti < t at which two or more discontinuities 
cross in the interior of fi. Observe that, since u'' is a front tracking solution constructed by 



the algorithm of Section 2.3 on the region [Ti,r] x [a,b], we can apply Lemma 3.3. Hence, it 



follows that the left and right states of the jumps occuring in u'^iji, •) satisfy condition (4 



because ( |3.20 ) guarantees that two (or more) adjacent rarefaction fronts of the same family cannot 
cross at time ti. We then solve backward the resulting Riemann problems applying again the 
above procedure. This determines the solution u'^{t, •) untill the time T2 < ti at which another 
intersection between its fronts takes place in the interior of Q, and so on (see figure la). 





figure la figure lb 

With this construction we define a front tracking solution u'^(t,x) on the whole 
region [(3/4)r, r] x [a, 6], that verifies the first equality in (4^), and corresponds to the boundary 
data = u^i-^a), = u''{-,b) e L°°([(3/4)T, r], T"). Clearly, the total number of wave-fronts 
in u'^{t, •) decreases, as t | (3/4)r, whenever a (backward) front crosses the boundary points 
X — a, X ~ b. Since (2.3) implies that the maximum time taken by fronts of u'^ to cross the 
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interval [a,b] is (6 — a)/A™'", the definition (4.4) of T guarantees that all the (backward) fronts 
of u"^ will hit the boundaries x = a, x = b within some time t' g](3/4)T, t[, which shows that 



also the second equality in (4.6) is verified, thus completing the proof. 



□ 



Lemma 4.2 Let T > be the constant in (4-4 )■ Then, for any pieeewise constant function 
jjv g L°°([a, 6], r''), and for every state uj e , there exists a front tracking solution 
u'^{t,x) of ( jj. jy on the region [0, (3/4) T] x [a,b], corresponding to some boundary data 
u'^, ul e L°°([0, (3/4)r],r''), so that 



(O, x) — u'^{x), 



u''{{3/4:)T, x) =LJ , yxe[a,b]. 



(4.10) 



Proof. The approximate solution u'^ is constructed as follows. By Remark 2.3 



for t e [0,T/4], we can define u^{t,x) as the restriction to the region [0,T/4] x [a, 6] of the 



front tracking solution to the Cauchy problem for (1.1), with initial data 

u{x) 



u^{a+) if X < a, 
vty {x) if a < a; < 6, 
u^ib-) if x>b, 



(constructed as in with the same type of algorithm described in Section p.3[ ). Observe 
that, since contains only fronts originated at the points of the segment {(0, x) ; x e [a,fo]}, 
because of ([2.5|), (4.4) these wave-fronts cross the whole interval [a, b] and exit from the boundaries 
X = a, X = b before time T/4 (see figure lb). Hence, there will be some state lo' G such that 



u-{T/A, x) 
Thus, introducing the intermediate state 



V a; e [a, b] 



(4.11) 



between lo' and w, we will define u'^{t,x), for t G [T/A, T/2], as the restriction to the region 



[T/4, T/2] X [a, b] of the approximate solution to the Riemann problem for (1.1), with initial data 

u''iT/4,x)^\ ^/ ' (4.12) 

I u(uj) if X > b , 

while, fortG [T/2, (3/4)T], we will let u''(t, x) be the restriction to the region [T/2, (3/4)T] x [a, 6] 



of the approximate solution to the Riemann problem for (1.1), with initial data 



i-{T/2,x) = 



u{uj) 



if 
if 



a; < a , 
X > a . 



(4.13) 



By the definition of w, and because of (2.3), (4.4), on [T/4, T/2] the solution of the Riemann 



problems with initial data ( 4.12 ) contains only wave- fronts originated at the point (T/4, 6), that 
cross the whole interval [a, b] and exit from the boundary x = a before time T/2. Similarly, still 
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by (|2J), for t G [T/2, (3/4)T] the solutfon of the Riemann problem with initial data ( [4.13D , 

contains only wave-fronts originated at (T/2,a), that cross the whole interval [a,b], and exit 
from the boundary x — b before time (3/4)T (see figure lb). Hence, u'^{t,x) is a front-tracking 
solution defined on the whole region [0, (3/4) T] x [a,b], that corresponds to the boundary data 
= «''(•, a), = ^{-,1) G L°°([0, (3/4)T],r'"), and verifies the conditions (p^). 



□ 



Step 2. Convergence of the approximate solutions. By Step 1, for a given (p G if (with p' 
as in ( 3.19| )), we have found a sequence of initial data u'^, and of boundary data mJ^, G U!^ , so 
that, letting u''{t, •) = El^.iu^ , u'^, u'j^) be the corresponding front tracking solution, there holds 



Li([a,6]). 



(4.14) 



By the same arguments used in the proof of Proposition we may assume that the initial 
and boundary data of each approximate solution u"^ have at most < v shocks for every 
characteristic family. Then, thanks to the Oleinik-type estimates ( 2.25| ), and because are 



uniformly bounded since they take values in the compact set ( 2.13 ), for every fixed £ > 0, there 
will be some constant > such that 



Tot.Var.{u''(t, • ) ; [a + £, b - e]] < 

rb-e 

\u'^{t, x) — u'^{s, x)\ dx < Ce\t 



G [e,T], 
s\ Vt, sG[£,r] 



(4.15) 



Hence, applying Helly's Theorem, we deduce that there exists a subsequence {u^^ }j>o that con- 
verges in L^([a,6], F) to some function Ue{t, •), for any t G [e, r]. Therefore, repeating the same 
construction in connection with a sequence e^, — > 0+, and using a diagonal procedure, we ob- 
tain a subsequence {u'^ {t, ■ )}i/'>o that converges in L^([a, b], T) to some function u{t, •), for any 
t G [0, r]. Then, by Theorem [2^, there holds ( |2.35| ), with Ua = u(-,a), ut = u{-,b) G W^, 



while (4.14) implies u(r, •) — ip, which shows (p G A(t). This completes the proof of Theorem 2.1 



□ 



We next establish the compactness of the attainable set ( ^.14 ) stated in Theorem 2.2. The 
proof is quite similar to that of Theorem 2.3]. We repeat it for completeness. 



Proof of Theorem 2.2 . Fix T > 0, and consider a sequence {u'^}i,>o of entropy weak solutions 



to the mixed problem for (LI) on CIt = [0, T] x [a, b] (according with Definition 2^), with a fixed 
initial data u G L°°([a, 5], F). Since all u"^ are uniformly bounded, and because of the Oleinik- 



type estimates (2.7)-(p.q), one can find, for every e > 0, some constant Ce > so that ( 4.15 ) 
holds. Thus, with the same arguments used in Step 2 of the previous proof, we can construct a 
subsequence {u" },y'>o so that, for any t G [0, T], u'^ {t, ■ ) converges in to some function u{t, •), 
which is continuous as a map from ]0, T] into L^([a,fe],F), and satisfies the entropy conditions 
( |2.?j )-( pTTo ) on the decay of positive waves. On the other hand, the weak traces '^'^ , '^'j^ of 
the fluxes f{u'^ ) at the boundaries x — a, x — b are uniformly bounded, and hence are weak* 
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relatively compact in L°°([0, T]). Thus, by possibly taking a further subsequence, we have 

K^'fa, in L°°([0,r]), (4.16) 

for some maps ^Pq, G L°°([0, T]). Notice that, by the properties of the Riemann invariants, 
the set /(r) is closed and convex, and hence also the weak limits 'i'a, take values in /(F). 



Moreover, since each u" is a distributional solution of (1.1)-(1.2) on fir, also the limit function 
M is a distributional solution of the Cauchy problem ( |l.l| )-(p^) on the region Ut- Then, setting 
Ua = f ° ^a, Ub = /^^ o v[/[,^ it follows that u is an entropy weak solution of the mixed 



problem (L1)-(L4) (with boundary data in U^) according with Definition p.2| , which shows that 



w(r, •) e A{T). This completes the proof of Theorem 2.2 



□ 



If we take in consideration only solutions to the mixed problem ( |l.l| )-( pr^ ) that are trajectories 
of the flow map E obtained in Theorem ^.3| (which, in particular, admit a strong trace at the 
boundaries x = a,x = b), we are lead to study the set of attainable profiles 



Ae{T) = {Et{u, 



Ub, 



7a, Kb GL°°([o,r],r)}, 



(4.17) 



Since AsiT) C A{T), and by the proof of Theorem 2.1, it clearly follows that the characterization 



of the set A{T) provided by the inclusions (|2.17|) - (|2.18D of Theorem ^ holds also for Ae{T). 
Concerning the compactness of the set Ae{T), observe that, given any sequence of exact solutions 
u'^{t, •) = Et{u'^ , u'^ , u^), V > 1, by Theorem p.3| one can find another sequence of approximate 



solutions v'^{t, ■) constructed by the front tracking algorithm of Section 2.2, so that 



Therefore, relying on the regularity property of a solution obtained as limit of front tracking ap- 



proximations provided by Theorem 2.5, with the same arguments used in the proof of Theorem 2.2 



one can establish also the compactness of the set Ae{T). 

Theorem 4.1 Under the same assumptions of Theorem 1, the set Ae{T) is a compact subset 
o/Li([a,6], r) for each T > 0. 
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